Lattice study on r] c2 and X(3872) 
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Properties of 2 h charmonium n C 2 are investigated in quenched lattice QCD. The mass of n C 2 is 
determined to be 3.80(3) GeV, which is close to the mass of D-wave charmonium ?/>(3770) and in 
agreement with quark model predictions. The transition width of nc2 — > jJ/ip is also obtained with 
a value of F = 3.8(9) keV. Since the possible 2 h assignment to X(3872) has not been ruled out by 
experiments, our results help to clarify the nature of X(3872). 

PACS numbers: 11.15.Ha, 12.38.Gc, 13.20.Gd, 14.40.Pq, 14.40.Rt 



I. INTRODUCTION 

Even though the charmoiumlike resonance A (3872) 
has been established for several years with Mj = 
3871.68 ± 0.17MeV and T x < 1.2MeV @, @, the very 
nature of it has not been fully understood till now. Its 
even C parity has been firmly established from its decay 
to J/ipp and to J/-07 [8]- Further analysis of its de- 
cay angular distribution also constrains its total quantum 
number J PC to be either 1 ++ or 2 h . The discovery of 
A (3872) has triggered quite a number of theoretical in- 
terpretations by assuming a quantum number such 
as the radial excitation of Xdi the DD* molecule [9l-tll|. a 
tetraquark state [l2l - fl4j ]. etc.; however, none of them can 
accommodate all the observed features of X(3872). The 
situation became more complicated when the BABAR 
Collaboration reported in 2009 that 2 h is more favored 
by the study of the decay angular distribution of the 
process X(3872) -» J/-07r + 7r~7r° [Hj]. In contrast, the 
same analysis by the Belle Collaboration claims that both 
the 1 ++ and 2 ^ assignments are consistent with their 
data 0, EH ■ Another controversial result comes from the 
measurements of the radiative decays of X(3872). For 
the decay mode X(3872) ->■ jJ/yb, the BABAR and Belle 
collaborations reported consistent measurements [l7l[l8| : 

Br(B ± V(3872) J ftT ± )Br(V(3872) J/ipj) 

= (2.8 ±0.8 ±0.1) x 10~ 6 (BABAR), 

Bt{B ± V(3872) J &T ± )Br(V(3872) J/ijr/) 

= (1.78t°; 4 *±0.12) x 10~ 6 (Belle). (1) 

With the world average value Br(B+ X(3872)K+) < 
3.2 x 10~ 4 , one can estimate the branch ratio 
Br(X(3872) J/ip'y) > 0.9% (BaBar) or 0.6% (Belle). 



* cheny@ihep.ac.cn 



However, for the decay mode X(3872) jip' , BABAR 
measured a 3.4 ±1.4 times larger branch ratio [T3|, but 
Belle found no evidence [llj]. This large discrepancy 
should be reconciled by further experimental measure- 
ments. 

Theoretically, if we are constrained to its charmonium 
assignments, Jf (3872) can be either the radial excitation 
of Xd (h sa Y: x'di or ^ ne 1 -^ ) 2 charmonium r) C 2 (if 

2 '"). The potential quark model predicts the mass of 
x'd to be 3925 MeV [19| , which deviates from the mass 
of X(3872) by about 50 MeV. There are also many lattice 
studies predicting a x' c i ma ss ranging from 3850 to 4060 
MeV [20 23], but with various uncertainties of their own, 
where the key difficulty is the challenging task of extract- 
ing the excited states. As for the r\ c i, the quark models 
usuallypredict the mass to be in the range 3770 to 3830 
MeV [191 . l24l . [25| . which is even further away from the 
mass of V(3872). This is also reinforced by recent lattice 
studies (and this work) . At any rate, the mass parameter 
should not be the unique criterion for the interpretation 
of X (3872); more information is definitely desired- for ex- 
ample, the radiative transition properties of x' e i an d ^2, 
which are theoretically accessible and hopefully can shed 
some light on the nature of X(3872). 

In this work, we will focus on the study of the prop- 
erties of rj C 2, such as its mass and radiative transition 
width to J I if). There are actually several phenomeno- 
logical studies on this topic |26l - |28| . but they are rather 
model dependent. In contrast, the lattice QCD approach, 
as a method from first principles, can provide information 
that is more model independent. An additional techni- 
cal advantage in the study of t? c2 on the lattice is that it 
is the ground state in the 2 ^ channel and is free from 
the uncertainty of the extraction of excited states. In 
view of the notorious bad signal-to-noise ratio for P and 
D wave states, we adopt the quenched approximation so 
as to obtain large enough statistics for precise physical 
quantities to be derived. As for the quenched approxi- 
mation, even though long-term experiences show that it 



2 



is safe for charm quark systems, and the resultant un- 
certainties can be small, we still take several steps to 
check this and be assured of our results. We first calcu- 
late the spectrum of the ground state charmonia, such 
as 15, IP states, and make sure that the experimental 
spectrum patterns are reproduced. As for the radiative 
transitions, we choose the transition mode of the tensor 
charmonium \ci to J/ip as a calibration of the systematic 
uncertainties of our formalism by comparing our result to 
the experimental value. After that, we continue to the 
study of the radiative transition of i] c2 to J/tp. All the 
lattice calculations are carried out on anisotropic lattices 
which are suitable to the study of heavy particles. The 
numerical techniques are standard: the mass spectra are 
extracted from two-point functions, and the multipole 
amplitudes contributing to the transition widths are de- 
rived from the calculation of relevant three-point func- 
tions with a local electromagnetic current insertion. We 
apply two anisotropic lattices with different lattice spac- 
ings to estimate the lattice artifacts owing to the finite 
lattice spacing. 

This work is organized as follows: The formalism for 
the calculation of radiative transition widths on the lat- 
tice is briefly introduced in Sec. HU In Sec. |TTT] are the 
numerical details, where the lattice setup, the extrac- 
tion of mass spectrum and transition form factors are 
explained, and the numerical results are presented. Sec- 
tion [IV] is the conclusion and discussion. The theoretical 
derivation of the multipole form factors is described in 
the Appendixes. 

II. FORMALISM 

As mentioned above, in this work we aim at the lat- 
tice calculation of the radiative transition rate of r\ c i to 
J/ip. The general radiative transition width of an initial 
particle i to a final particle / is 

X l M n,r„r 7 | 2 , (2) 

where q — pi — pj is the decay momentum with the mass- 
on-shell value \q\ = (M 2 - M|)/(2M l ), M, and M f are 
the masses of the particles i and /, and M r ^ r , r is the 
transition amplitude with ri,r/,r 7 being the polariza- 
tions of i, /, and the photon, respectively. To the lowest 
order of QED, the amplitude M is expressed explicitly 
as 

M rtm =^(^r 7 )</(? / ,r / )|^(0)|i(ft,r i )), (3) 

where e*(<f, r 7 ) is the polarization vector of the pho- 
ton, and (f(p', ^/)b'^ m (0)|i(p, r^)) gives the on-shell ma- 
trix elements of the electromagnetic current je m (x) = 
i/jQ^f ,i ip(x) between the i and / states. [Here ijj refers 



to an array of all the contributing quark flavors, such as 
u, d, s, c, . . ., and Q is a diagonal matrix of quark elec- 
tric charges, say, diag(Q) = Q u ,Qd,Qs,Qc, ■ ■ ■■] The 
hadronic matrix element can be derived directly from the 
lattice QCD calculation of the related three-point func- 
tions: 

T^t(p f ,q, t, t>) = e-^e-^ x 

(Ottg,t)jj^(MO$@,0)), (4) 

where O^/n are the interpolating fields for the particles i 
and /, with the indices to, n referring to different spatial 
components for spin nonzero states. The explicit deriva- 
tion can be expressed as 

r^(p f ,q,t,t') = £ e -%*e-(^/)*' 
n,r f 

x zL(j>f, r f)Zn(pun) 

2E t 2E f 

(t',t-t' ->■ oo), (5) 

where Zfyf are the matrix elements like Z*(px,rx) = 
(0|O^|X(px, fx)), which can be derived from the rele- 
vant two-point functions, 

r% mn @x,t) = £e'*-*<Q£(*,i)0*t(ff,o)) 

—x 




J2(0\O*\X(p x ,rx)) 



rx 

x (X(p x ,r x )\O^\0) (t->oo). (6) 

On the other hand, in the Minkowski space-time, the 
matrix elements (/(p/, Tf)\j^ m (0)\i(pi, Tj)) can be gen- 
erally expressed by several Lorentz-invariant form fac- 
tors Fk(Q 2 ) and Lorentz-covariant kinematic factors 
a k(j>i,Pf) through the multipole decomposition, 

(f(p f ,r f )\]Um^r l ))=Y<<^P^ F ^ 2 ^ ( 7 ) 

k 

where pij are now the four-momenta of particles i and 
/, and Q 2 is the squared transfer momentum Q 2 — 
— (pi —pf) 2 . Obviously, the concrete form factors Fk(Q 2 ) 
and the explicit expressions of the kinematic factors ctk 
depend on the properties of the particles i and /, and 
therefore should be worked out case by case. Finally, the 
decay width with an on- shell photon (Q 2 = 0) can be 
expressed as 

r(^ 7 /)«^F 2 (0). (8) 

fe 

So the key problem in this work is to reliably extract 
these form factors through the lattice calculation of the 
relevant hadronic two-point functions and three-point 
functions described above. 
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TABLE I: Relevant input parameters for this work. The spa- 
tial lattice spacing a a is determined from r^ 1 = 410(20) MeV 
by calculating the static potential. 



p 




a s (fm) 


ia s (fm) 


L 3 x T 


N con f 


2.4 


5 


0.222 


1.78 


8 a x 96 


1000 


2.8 


5 


0.138 


1.66 


12 3 x 144 


1000 



III. NUMERICAL DETAILS 

We use the quenched approximation in this study. The 
gauge configurations are generated by the tadpole im- 
proved gauge action [2!| on anisotropic lattices with the 
temporal lattice much finer than the spatial lattice, say, 
£ = a s /a t 1, where a s and at are the spatial and 
temporal lattice spacings, respectively. The much finer 
lattice in the temporal direction yields a higher resolution 
to hadron correlation functions, such that the masses of 
heavy particles can be tackled on relatively coarse lat- 
tices. We have two anisotropic lattices (L 3 x T — 8 3 x 96 
and 12 3 x 144) with £ = 5. The relevant input pa- 
rameters are listed in Table HI where the lattice spac- 
ings, say, a s = 0.222(2) fm for the coarser lattice and 
a s = 0.138(1) fm for the finer lattice, are determined 
from r^ 1 = 410(20) MeV by calculating the static poten- 
tial. For each lattice, we generate 1000 configurations, 
each of which is separated by 500 heat-bath updating 
sweeps to avoid the autocorrelation. For fermions, we 
use the tadpole improved clover action for anisotropic 
lattices [3(J. The parameters in the action are tuned 
carefully by requiring that the physical dispersion rela- 
tions of vector and pseudoscalar mesons are correctly re- 
produced at each bare quark mass [3l[. The bare charm 
quark masses for the two lattices are set by the physical 
mass of J/tp mj/^ = 3.097 GeV. 

In this work, we only consider the connected diagrams 
in the calculation of two-point and three- point func- 
tions. The contribution of the disconnected diagrams 
is assumed to be small for charmonium states due to the 
OZI suppression. 

A. Ground-state charmonium spectrum 

As the first step, we carry out a careful study on 
the ground-state charmonium spectrum, which can illus- 
trate to some extent the systematic uncertainties due to 
the quenched approximation. For the states ry c (0 h ), 
JM1 — ), h c (l+~), Xco(0++), and X d(l + +), we adopt 
the conventional quark bilinear operators like cFc, with 
r = 75, 7i, (Tij, 1, and 75 7.;, respectively. For the tensor 
mesons x c2 (2 ++ ) and 77^2 (2 h ) , since there are not quark 
bilinear operators, we build the corresponding operators 
by combining the quark bilinear operator with either the 
spatial gauge-covariant derivatives Z3j or the color mag- 
netic field strength operator Bi , which is built from Wil- 
son loops. It is known that the spin J = 2 states in the 



continuum correspond to both the T2 and E irreducible 
representations (irreps) of the cubic point group O on 
finite lattices, so the interpolating field operators of the 
two irreps are constructed for the tensor charmonia. For 
example, the To operator for the y C 2(2 ++ ) state is taken 

^ y ^ y / y 

as \tijk\c'jj D kC where D = D — D , and the E operator 
is also built to check the restoration of the continuum 
rotation symmetry. 

We will emphasize the choice of the operators for the 
2 h state, which is the major object of this work. The 
situation for the ry C 2 meson is a little bit more compli- 
cated. We try first three types of operators, such as 

^ y 

\e ljk \c{x)Y,j D k c{x) (F>-type), 

^ y ^ y 

\eij k \c(x)j 5 D j D k c(x) (£>D-type), 
\eij k \c(x)jjB k c(x) (F-type), 

where only the T2 operators are presented (F operators 
can be built similarly, and the details can be found in 
Ref. H). 

It is known that the signal-to-noise ratios of the cor- 
relation functions are always bad for P-wave and D- 
wave states. To circumvent this difficulty, we adopt the 
Coulomb-gauge fixed wall source techniques in the calcu- 
lation of the spectrum. The configurations are fixed to 
the Coulomb gauge first, then the charm quark propaga- 
tors are calculated with uniform wall source vectors. For 
the spin J = 0, 1 states, the point-sink wall source cor- 
relation functions can be constructed straightforwardly 
with these propagators. For the tensors, we use the 
F-type operators as the wall source, which means that 
additional inversions should be carried out with wall 
sources multiplied by the local color field strength op- 
erators Bi{x). On the other hand, since the gauge is 

fixed, the gauge-covariant derivative operator D is re- 
placed by the direct derivative operator V = V - ^ in 
the practical calculation. 

The masses of IS and IP charmonium states can be 
neatly derived with the standard data analysis, however, 
the situation for the 2 h channel is very strange. Fig- 
ure [1] shows the effective masses of various correlation 
functions of this channel at j3 = 2.8. It is seen that 
the effective mass of the F-type point sink and -F-type 
wall source correlator (F — F) saturate at a plateau with 
the best-fit mass 4.43(8) GeV, while that of the DD-type 
point sink and F-type wall source correlator (DD — F) 
goes lower and does not show a perfect plateau. Intu- 
itively, a mass of 4.4 GeV is too large for the 2 ^ ground 
state charnonium. Thus what one can infer from these 
behaviors is that the F-type operator couples predom- 
inantly to a higher state but little to the conventional 
charmonium; in the mean time, there must be a lower 
state which can be accessed by the DD-type operator 
but whose spectral weight is relatively small due to the 
F-type wall source. To check this and to dig out the de- 
sired 2 ^ charmonium state, we try instead another wall 
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F(wall)-F(point) ' — ' — ' 

F{wall)-DD{point) ' — « — ' 

6-5 - Q(wall)-F point ' — ' — ' ' 

t Q(wall)-DD{point) ' — « — • 

fi * 4.43(8)GeV 

3.79(3)GeV 

5-5 ■ + * 4 T - 

5 - "% t 



5 10 15 20 25 30 35 40 
t/a 

FIG. 1: The 2" + effective masses of F - F, DD - F, F - Q 
and DD — Q correlators at /3 = 2.8 are plotted for illustration. 
F—F and F—Q effective masses lie on each other and saturate 
to the same plateau with the best-fit mass M = 4.43(8) GeV 
(indicated by the red line with the jackknife error band). The 
DD — Q effective mass has a plateau at M — 3.79(3) ( blue 
line with the jackknife error band ). In contrast, DD — F 
effective mass does not show perfect plateau, but it evolves 
gradually from the upper plateau to the lower. 



source operator (T2 irreps for example), 

|e tffc | c a (x,0h 3 c' a (x,0)c> b (y,0)X k c b (z,0) (Q-type), 

where = eijfc<7y and d stands for a quark field with 
the same mass as that of the charm quark but a different 
flavor. For simplicity, we call this operator Q-type in the 
context. With this type of wall source operator, the effec- 
tive masses of the F-type point sink correlator (F — Q) 
and the DD type point sink correlator (DD — Q) are 
also plotted in Fig. [TJ where one can find that the mass 
plateau of the F — Q correlator coincides with that of 
the F — F correlator within errors, while the effective 
mass of the DD — Q correlator shows a very nice plateau 
with the best-fit mass 3.79(3) GeV. Since the lower state 
has a mass close to the potential model prediction of 
2 *~ charmonium and the higher state is much heavier, 
we assign the lower state to the conventional 1 £>2 char- 
monium state rj C 2- This assignment can be reinforced 
by comparison with the established l 3 _Di charmonium 
state ■0(3770): They are both D-wave charmonia and 
are therefore close in mass; the small mass splitting can 
be attributed to the different spin-spin and spin-orbital 
interactions. 

The whole spectrum of the lowest-lying charmonium 
states we extracted in this work is illustrated in Fig. [2] 
and listed in Table HT1 where the experimental values are 
also given for comparison. Since we have only two lattice 
spacings, we would not carry out a serious extrapolation 
to the continuum limit, but we show all the results, from 
which one can see that the effects of the finite lattice 
artifacts and the quenched approximation are not that 
important. 



TABLE II: Listed here are the masses of the lowest-lying char- 
moinum states extracted from the two lattices (/3 = 2.4 and 
/3 — 2.8) in this work. The experimental results 5] and the 
nonrelativistic quark model predictions [lj| are also given for 



comparison. 


meson 


j pc 


M(2.4) 


M(2.8) 


Expt. 


QM 


VcCSo) 
J/V>( 3 Si) 


(r 

r 


-+ 


2.989(2) 
3.094(3) 


3.007(3) 
3.094(3) 


2.981 
3.097 


2.982 
3.090 


hoCPo) 
Xco( 3 Po) 

XclfPl) 
Xc2( 3 P 2 ) 


1+- 

0++ 

1++ 

2++ 


3.530(35) 
3.472(34) 
3.508(50) 
3.552(17) 


3.513(14) 
3.431(30) 
3.499(25) 
3.520(15) 


3.526 
3.415 
3.511 
3.556 


3.516 
3.424 
3.505 
3.556 


r( 3 D 1 ) 
vc2Cd 2 ) 


r 

2" 


-+ 


3.777(30) 


3.789(28) 


3.770 


3.785 
3.799 





P=2.4 1 1 




P=2.8 1 1 




Exp ' — x — 1 




V (3770) D 




XcO X C 1 Xc2 h c ^2 






JA|< 




_% -H- 









1 S 3 S 1 3 P 3 P, 3 P 2 ^ 3 D, 1 D 2 



FIG. 2: IS, IP, and ID charmonium spectrum. The red 
boxes illustrate the results for /3 = 2.4, and the blue ones for 
ft = 2.8. The experimental value are also plotted with points 
for comparison. 



The goal of the spectroscopy study in this work is 
twofold. First, the physical spectrum of experimentally 
established charmonium states can be well reproduced in 
our formalism. This gives us confidence in our predici- 
ton of the r\ C 2 mass. Second, the practical study finds 
that the DD-type operator is preferable for producing 
the 2 charmonium. Therefore, in the study of its ra- 
diative transition, we choose the DD-type operator for 
r/ C 2 in the calculation of the related three-point functions. 



B. Renormalization of the vector current 

In the quenched approximation, since there are no 
sea quarks, the electromagnetic current contributing to 
the radiative transitions of charmonia involves only the 
charm quark, say, j cm {x) = Qd^ix) with j^(x) — 
cq^c(x), which is the one we adopt in this study. It is a 
conserved vector current and need not be renormalized in 
the continuum. However, on a finite lattice, it is not con- 
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and Z { y ] of 



TABLE III: The renormalization constants Z 
the spatial and temporal components of the vector current for 
P — 2.4 and ft — 2.8 lattices. Two momentum modes, (0,0,0) 
and (1,0,0), are used for the derivation. 



p 


zf (0,0,0) 


Zf (1,0,0) 


4 S) (i.o,o) 


2.4 
2.8 


1.288(5) 
1.155(3) 


1.299(11) 
1.159(3) 


1.388(15) 
1.110(7) 



served anymore due to the lattice artifact and receives a 
multiplicative renormalization factor Zy(a s ). Following 
the scheme proposed by Ref. [33], Zy(a s ) is extracted us- 
ing the ratio of the r\ c two-point function and the related 
three-point function evaluated at Q 2 = 0, 



e(j>) m^(p,p,f,t) ' 



where the factor 1 /2 accounts for the effect of the tempo- 
ral periodic boundary condition, and the superscript /i of 
Zyicis) is used to differentiate the temporal component 
from the spatial ones, since they are not necessarily the 
same due to the anisotropic lattices we use. Figure IIII Bl 
plots Zy(t) with respect to t for the two lattices. Zy's 
are extracted from the plateaus and the values are listed 
in Table IIII1 Obviously, the spatial components Zy ^ (a) 
deviate from the temporal ones by a few percent. This 
deviation can be attributed to the imperfect tuning of 
the bare velocity in the fermion action. In this work, 

(s) 

only Z v s enter the calculation since only the spatial 
components of the vector current are involved in the ex- 
traction of the form factors. 



C. Three-point functions and form factors 

With the prescriptions discussed above, we now give 
a brief description of the calculation of the three-point 
functions. In practice, we use local sink and source op- 
erators for the initial and the final states, and insert the 
vector current j^ix) = 07^0(2; ) only on the quark line. 
(The current insertion on the antiquark line is numeri- 
cally equivalent and is taken into consideration by mul- 
tiplying by a factor of 2 in the final result.) The three- 
point functions contributed by the connected diagrams 
(disconnected diagrams are neglected) are calculated by 
using the standard sequential source technique. (One 
can refer to Refs.[33l. l34j for the details.) In order to in- 
crease the statistics, we repeat the same calculations T 
times (where T is the temporal lattice size) by setting 
a point source on a different time slice each time. With 
the related two-point functions calculated accordingly, a 
straightforward way to extract the interested matrix el- 
ements (f(p', r/)|j M (0)|i(p, r,)) is to fit the three-point 
function and two-point function simultaneously accord- 
ing to Eq. ([5]) and Eq. ©. However, it is known that the 



1.6 
1.55 

1.5 
1.45 

1.4 
1.35 

1.3 
1.25 

1.2 



1.25 



Z,(0,0,0) 
Z t (1,0,0) 
Z,(1,0,0) 



1.05 



10 15 20 25 30 35 40 45 

t/a 



2^(0,0,0') 



(1,0,0) 



ilUJ 



20 



30 40 

t/a 



50 



70 



FIG. 3: The renormalization constant Zy of the vector cur- 
rent. The upper panel is for /3 = 2.4 and the lower one for 
/3 = 2.8. The dots are the simulation data with jackknife 
errors, and the lines show the fit results. 



excited states contribute much to two-point and three- 
point functions when the time ranges t and tf — t are 
not large enough. This situation is more serious for local 
operators, so it is not trivial to isolate the contribution of 
ground states. A way around this is to employ the ratios 
of correlation functions, which can suppress the contri- 
bution of excited states substantially. For this purpose, 
we introduce the functions 



Rf>(t) = rW(pf,q,t f ,t) X 



2E l Tf\p l ,t f -t) 



\Tf\pXt)Tf\p ll t f )\vf{p- h t f 



2E f rf' 



(9) 



which should be insensitive to the variation of t in 
a time window, so that the desired matrix element 
(/(P/j r /)|j M (0)|V'N(Pi: r i)) can be extracted from the 
plateau. 

In the data analysis, we divide the 1000 configurations 
into 100 bins and use each bin average as an indepen- 
dent measurement. For the resultant 100 bins, we use 
the one-eliminating jackknife method. Since the energies 
Eij can be determined very precisely from the two-point 
functions, they are treated as known parameters in the 
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above equation. Practically, is fitted by the func- 

tion 

fl"(t) = (f(Pf,r f )\fiO)\^m,r t ))+Sf(t) (fO) 

where the additional term Sf(t) — aer bmt accounts for 
the residual contribution of excited states. Thus we can 
obtain a jackknife ensemble of the matrix elements. The 
second step of data analysis is to extract the form factors 
that enter the calculation of decay widths. Since these 
matrix elements can be expressed in terms of form factors 
through the multipole decomposition, 

(f(p f ,r f )\j»(0)m,ri)) =J2<(P^Pf) P k(Q 2 ), (11) 

k 

and a^(j>i,pf) are theoretically known kinematic func- 
tions, the form factors Fk(Q 2 ) can then be derived 
straightforwardly. Taking into consideration the contri- 
bution of the current insertion on the antiquark line, the 
electric charge of the charm quark Q c — 2/3, and the 
renormalization constant of the spatial components of 
the current operator Zy , Fk is related to Fk of Eq. ([7} 
as 

F k {Q 2 ) = 2x^ex Z ( v s) F k (Q 2 ). (12) 
With this in mind, in the following context, we omit the 

A (s) 

hat of F and insert Z v implicitly in possible expressions. 

In order to take good care of the correlation between 
the form factors, we carry out correlated minimal % 2 fits 
with the jackknife covariance matrix built from the jack- 
knife ensemble of the matrix elements. On the other 
hand, for a specific Q 2 , there may be several symmetric 
copies of the matrix elements with the same value of . 
These copies are averaged over to increase statistics. 

In the following subsections, we present first the cal- 
culation of the process Xc2 — > jJ/ip to see how precisely 
the form factors- and thereby the transition width- can 
be derived, and then the results of r\ c2 — > ^ip. 

D- Xc2 — ► yJ/ii transition 

The Minkowski space-time matrix elements for this 
transition can be expressed in terms of form factors as 
follows: 

(V(P v ,XvW(0)\T(Pt,Xt)) = <i?i(Q 2 ) 
+a%M 2 {Q 2 ) + a%E 3 {Q 2 ) + <d(Q 2 ) + a»C 2 (Q 2 ) 

(13) 

where V stands for the 1 vector meson J/tp, T stands 
for the 2 ++ tensor Xc2, and a? are Lorentz covariant kine- 
matic functions of pv and pr (and specific polarizations 
of V and T), whose explicit expressions are tedious and 
omitted here. Although a J = 2 representation of the 
rotational symmetry in the continuum breaks into the 
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FIG. 4: The extracted form factors E 1 (Q 2 ), M 2 (Q 2 ), and 
Es(Q 2 ) are plotted versus Q 2 for the two lattices of (3 = 2.4 
(the upper panel) and f3 = 2.8 (the lower panel), respectively, 
where the points are the simulation data, the line the fit func- 
tion, and the error bands the jackknife ones. The PDG values 
of Ei (0) and M2 (0) are also plotted for comparison. 



E and T 2 irreducible representations (irreps) of the lat- 
tice spatial symmetry group O, we find that this break- 
ing effect is small in our work, as is manifested by the 
near degeneracy of the masses and spectral weights of 
the ground states in these two irreps when we study the 
relevant two-point functions. Thus, we assume that the 
rotation symmetry breaking is also negligible for the re- 
lated matrix elements, and we carry out the multipole 
decomposition on the basis of E ® T 2 , which is equiva- 
lent to the J — 2 basis up to an orthogonal transforma- 
tion. One can find the detailed decomposition procedure 
in Appendix [3] and may also refer to Refs. [33|, |35|. In 
the practical study, we set T to be at rest and let V move 
with different spatial momenta p — 2nn/L. The 27 mo- 
mentum modes of ft = (^1,712,^3) ranging from (0,0,0) 
to (2, 2, 2) are calculated for V. 

The transition width of \ c2 — > 7J/V' for an on-shcll 
photon (Q 2 = 0) involves only the form factors £a(0), 
M 2 (0) and E 3 (0), or explicitly, 

T( Xc2 -> 7 J/tf) = ^^(| Sl (0)| 2 + |M 2 (0)| 2 + |i?3(0)| 2 ) 

Xc2 

(14) 
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where \q\ = {M 2 a2 - Mj,^)/2M Xc2 is the decaying en- 
ergy of the photon, and a — e 2 /(An) is the fine struc- 
ture constant. Since our simulation data are obtained at 
Q 2 0, these on-shell form factors should be interpolated 
to Q 2 = 0. In doing this, we adopt the fitting functional 
form inspired by the nonrelativistic quark model 33] , 



F k (Q 2 ) = F k (0)(l + X k Q 2 )e 



(15) 



which has been applied successfully in previous works. 
Here F k (0), X k , and ft are the parameters to be fitted 
though a correlated \ 2 fitting procedure where the co- 
variance matrix is constructed using the one-eliminating 
jackknife method. Plotted in Fig. @] are the extracted 
form factors Ei(Q 2 ), M 2 (Q 2 ), and E 3 (Q 2 ) versus Q 2 for 
the two lattices of f3 = 2.4 (the upper panel) and /3 = 2.8 
(the lower panel). The data points are the simulation re- 
sults, and the lines are the fit function with the jackknife 
error bands. One can find that the data are very precise 
owing to the high statistics, and the fit errors are also 
very small. We also carry out a simple polynomial fit with 
respect to Q 2 , F k {Q 2 ) = F k (0) + aQ 2 + c 2 Q 4 , and get 
consistent results within errors. Table HVl lists the results 
of the interpolation, where the continuum limit extrapo- 
lation is also given. It is seen that the electric dipole (Ei) 
contribution is dominant in the transition \ c2 lJ/^P, 
while the contribution of the magnetic quadrupole (M 2 ) 
is drastically suppressed, as depicted by the ratio 



a 2 



M 2 (0) 



^i(0) 2 + Af 2 (0) 2 +£;3(0) 2 



(16) 



for which we get a result a 2 = —0.107(3) for /3 = 2.4 
and a 2 = —0.082(7) for j3 = 2.8. After a linear extrap- 
olation in a 2 , we get the value in the continuum limit 
a 2 = —0.067(7), which is consistent with the PDG data, 
where a 2 = -0.100 ± 0.015 ||. The contribution of the 
electric octupole E 3 is far smaller. For the ratio 



03 



£3(0) 



y/E^oy + AUoy + E^oy 



(17) 



we obtain a 3 = 0.007(2) for /3 = 2.4, a 3 = 0.003(4) for 
/3 = 2.8, and the continuum limit a 3 = —0.003(6), which 
are also compatible with the PDG data a 3 = 0.016 ± 
0.013 0- If we focus on E%, we get the fitting parameters 
ft and Ai, 



ft = 0.431(5) GeV 
\i = -0.285(2) GeV -2 



for (3 = 2.4 and 



ft = 0.395(4) GeV 
Ai = -0.336(3) GeV" 2 



(18) 



(19) 



for (3 = 2.8. 

Using the interpolated form factors F k (0) and taking 
the fine structure constant a — 1/137, the transition 



TABLE IV: Listed here are the results of the interpolated 
form factors -Ei(O), M2(0), and £3(0), as well as the transition 
widths. The continuum limits are also given. All the results 
are in physical units. The widths can be compared with the 
PDG data F = 380(30) keV @]. 



p 


Si (GeV) 


M 2 (GeV) 


S 3 (GeV) 


T(keV) 


2.4 


2.04(2) 


-0.218(4) 


0.014(3) 


347 ± 20 


2.8 


2.08(2) 


-0.171(10) 


0.005(8) 


352 ± 11 


Cont. 


2.11(2) 


-0.141(15) 


-0.007(12) 


361 ±9 



width can be calculated directly. As shown in Table IIV1 
the partial decay width of T(\c2 -> ^ 7) is predicted 
to be 347 ± 20 keV or 352 ± 11 keV for the two lat- 
tices, respectively. The continuum extrapolation gives 
r = 361 ± 9 keV. All these results can be compared with 
the PDG average of 380(30) keV. The agreement with 
experimental data of the Xc2 J H- 3 7 transition indi- 
cates that our method for the \c2 — > J/ip transition is 
reliable. Then we can turn to transition the rj c2 — > "fj/t/;. 



E. 



7/^7 transition 



The general Lorentz decomposition of the Minkowski 
matrix elements responsible for the transition rj c2 — > 
7 <//?/> can be expressed as 

(V(pv,\vW(Q 2 )\T(pt,\t)) =a(Q 2 )A» 
+b{Q 2 )B^ + c{Q 2 )a i + d(Q 2 )D' 1 + e{Q 2 )E»{2Q) 

where T stands now for the tensor meson r] c2l a(Q 2 ), 
b(Q 2 ), c(Q 2 ), d(Q 2 ), and e(Q 2 ) are Lorentz-invariant 
scalar functions of Q 2 , and A^, B^, C M , are kine- 
matic functions whose explicit expressions can be found 
in Appendix [Bj With the multipole decomposition, the 
matrix elements can be also expressed in terms of form 
factors Mi, E 2 , M 3 , and C 2 : 



(V(pv,XvW(0)\T(p T ,X T )) 
+ia% E 2 (Q 2 )+ia^ M 3 (Q 2 ) 



ia%C 2 (Q 2 ) (21) 



where are also kinematic functions which can be ex- 
pressed in terms of the kinematic functions in Eq. (|20l) . 
(See Appendix [Bj) With real photons in the transition 
Vc2 jJ/^j only three multipoles are contributing: the 
magnetic dipole (Mi), the electric quadrupole (E 2 ), and 
M3. The transition width is written as 



16a\q\ 

45M 2 c2 



(|M!(0)| 2 + |£; 2 (0)| 2 +|A/3(0)| S 



(22) 

Since they are calculated at Q 2 ^ 0, the multipole am- 
plitudes should be interpolated to Q 2 — 0. The form fac- 
tor C 2 (Q 2 ) corresponds to the emission of longitudinal 
photons and does not contribute at Q 2 — 0. In extract- 
ing the amplitudes, we take the standard procedure as 
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described in Sec. II. The three-point functions are calcu- 
lated by setting the tensor at rest and making the vector 
moving. In analogy with the 2 ++ case, the effect of rota- 
tional symmetry breaking between T 2 irreps and E irreps 
is found to be small in this case and is neglected in the 
data analysis. The form factors Mi(Q 2 ), E2(Q 2 ), and 
Af 3 (Q 2 ) with various Q 2 are extracted jointly by a corre- 
lated fitting with a one-eliminating jackknife covariance 
matrix, and the results are illustrated in Fig. [5] as data 
points with jackknife errors. In the following discussion, 
we will focus on the interpolation procedure. What is 



Mi (ft) 



£ 2 (ft) 



M 3 (ft) 



It is seen that each multipole form factor can be expressed 
as a series of ^./(mym 2 ^) with a prefactor \/fi/rriT. In 
the rest frame of T (as is the case in our calculation), the 
expression of ft is simplified as ft = (m,T\pv\) 2 , such that 
ft / '(rriym 2 ^) = v 2 , with v = \pv\/mv being the spatial 
velocity of V. The convergence of the series in v is guar- 
anteed if the form factors a, 6, c, d, e are not singular in 
Q 2 , since v < 1 (For our calculation, the largest value of 
v is approximately 0.5). So, for the decaying T at rest, 
we have the simplified expression of the form factors M\ , 
E 2 , and M 3 : 

Mi = \p v \(A 1 (Q 2 ) + B 1 (Q 2 )v 2 + C 1 (Q 2 )v 4 + 0(v 6 )) 
E 2 = \pv\(A 2 (Q 2 ) + B 2 (Q 2 )v 2 + C 2 (Q 2 )v 4 + 0(v 6 )) 
M 3 = \pv\(B 3 (Q 2 )v 2 +C 3 (Q 2 )v 4 + 0(v 6 )). (25) 

With these expressions, the following information can be 
inferred: (i) The desired \pv\ prefactor accounting for 
the P-wave decay of r\ c2 — > jJ/ip is explicitly derived, 
(ii) The leading contribution to Mi and E 2 is of order 
0(1), while that of M 3 is of order 0(v 2 ). For the case 
of this study, since f max ~ 0.5, it is reasonable that the 
nonsingular A4 , Bi and C\ can be expanded with respect 



interesting is the relation between the two sets of form 
factors. We should first mention that the form factors, 
given the fact that they are functions of Q 2 , can also be 
written in terms of another Lorentz invariant variable, ft: 

ft = (pv ■ Pt) 2 — m v m T 

= h(mv + m T f + Q 2 ][(my - m T f + Q 2 }. (23) 

Thus the two sets of the form factors are related to each 
other as follows, 



I 

to v, such that we can take the following functions to do 
the interpolation: 

F t {v) = Av + Bv 3 + Cv 5 + 0{v 6 )(F % -> Mi, E 2 ) 
Fi(v) = Bv 3 + Cv 5 + Dv 7 + 0(v 9 )(F t -> M 3 ), (26) 

and the on-shell amplitudes M\(Q 2 = 0), E 2 (Q 2 — 
0), and M 3 (Q 2 = 0) can be reached by Fi(vo) with 
vq = [m3p — m v ) / (2m,Tmv)- The extracted form fac- 
tor and the interpolation are shown in Fig. [SJ where 
the data points are the simulated results with jackknife 
errors. One can see that at v = [corresponding to 
Q 2 = —{uiT — mv) 2 ~ 0.5 GeV 2 ] the form factors Mi, 
E 2l and M 3 are surely consistent with zero. The fits 
using Eq. (|26[) are also shown as curves with jackknife 
error bands. The interpolated values of these form fac- 
tors at Q 2 = for both /3 = 2.4 and /3 = 2.8 are listed 
in Table |Vj where the resultant transition widths and 
the corresponding continuum limits are also given. It is 
surprising that, for both lattices, the obtained |M 3 | is 
unexpectedly large and comparable to M\. This may be 
qualitatively attributed to recoiling effects of the charm 
quark or charm antiquark by emitting the hard photon 
with an energy £L ~ 0.6 GeV in this transition, which 



.Vft 

rriT 



(a my + a m>T ~ 2c tjit) + 



2a my — 3c mx — 4d mhrnx + 6e mym\ 



4Vl5 



ft 



-2a my + 3c tot 



16V15 



ft 



O 



ft 



2 2 
mfpTYly 



.Vft 

tot 



ft 



/ 3 . 2a mv — c tot — 4c? mhrnx + 2e mvm' 

- (a tot — a mv) H 7= \ 

4 v v ' 4^3 V mlm? 



+ 



-2a mv + c tot 



16>/3 



ft 



2 2 

TO T TOy 



o 



ft 



l/rp 1 1 ' J ~y 



.Vft 

tot 
f 01 



-a mv — c tot + 2d m\mx + 2e mvm 



ft 



<ji~p 1 1 by 



ft 



a mv + c Tot 
4VT5 



ft 



Vrpllly 



(24) 
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may result in large form factors d(Q 2 ) and e(Q 2 ) (see the 
discussion below) . In contrast to the mild dependence of 
Mi and M3 on the lattice spacing, the form factor E2 
is very sensitive to the lattice spacing. The reason for 
this is unclear and under investigation. Anyway, after a 
naive continuum extrapolation using the data from the 
two lattices in this work, we get the continuum results of 
the form factors as follows: 



Mi = 0.104(10)GeV, 
E 2 = -0.071(20)GeV, 
M 3 = -0.132(10)GeV. 



(27) 



Applying these results to Eq. (f2"2"j) . the transition width 
of r/ C 2 — > jJ/tp is predicted to be 



T(rj c2 -» jJ/ip) = 3.8 ± 0.9 keV. 



(28) 



There have also been several phenomenological studies 



«• -0.5 




FIG. 5: The r/ C 2 — J/ip transition form factors Mi(Q 2 ), 
E 2 (Q 2 ), and M 3 {Q 2 ) are plotted versus Q 2 for the two lat- 
tices of f) = 2.4 (the upper panel) and /3 = 2.8 (the lower 
panel), respectively. The points are the simulation data, and 
the lines illustrate the fit function with jackknife error bands. 



on this transition, one of which is in the framework of the 
light- front quark model [27| , where the on-shell transition 



TABLE V: Listed here are the interpolated values of the form 
factors Mi, E 2 , and M 3 at Q 2 = for both j3 = 2.4 and /3 = 
2.8. The resultant transition widths and the corresponding 
continuum limits are also given. 



p 


Mi(GeV) 


£ 2 (GeV) 


M 3 (GeV) 


T(keV) 


2.4 


0.133(13) 


0.111(17) 


-0.093(9) 


4.4 ±0.9 


2.8 


0.115(11) 


-0.0007(14) 


-0.117(9) 


3.1 ±0.6 


Cont. 


0.104(10) 


-0.071(20) 


-0.132(10) 


3.8 ±0.9 



amplitude is decomposited as 

(V(pvAvW(0)\V2(pt,\t)) = 
~1h^ pa V T v V V p £{VT, X T )4(f v , Ay) 
+(/ 2 + f 3 )e^p^eUpv, \v)e a0 (p T , \t)pZp% 
+2he^plp v p ei(p T , \ T )p^e* a (p v , X v )pl] , (29) 

and the effective couplings are determined to be 



/1 = -0.0140(2)GeV _1 , 

f 2 = 0.146(3)GeV~ 3 , 

h = -0.092(l)GeV" 3 , 

U = 0.0180(l)GeV" 3 . 



(30) 



Since this decomposition is equivalent to Eq. (|2"0"|) by the 
relation 

c(0) = 2/ 1; d(0) = -(/ 2 ± / 3 ), e(0) = -2/4 (31) 

[It should be notified that a(Q 2 ) and b(Q 2 ) are equal 
to zero when Q 2 = because they are proportion to 
C 2 (Q 2 )-], the corresponding multipole amplitudes can be 
calculated from Eq. (j2"4")l as 

Mi = 0.079(2)GeV, 

E 2 = -0.086(2)GeV, 

M 3 = -0.125(3)GeV, (32) 

which gives a width of r=3. 54(12) keV. Taking into con- 
sideration the uncertainty of the choice of parameters 
such as the charm quark mass m c and the wave function 
parameter, etc., one can find that the lattice results and 
the LFQM results are surprisingly in excellent agreement. 
On the other hand, with the values in Eq. (|3H)) . we find 
that the coefficients Bi(Q 2 ) of the v 2 term in Eq. ([23)) are 
surely much larger than Ai(Q 2 ) at Q 2 = so as to com- 
pensate for the suppression of v 2 . This explains to some 
extent the fact that the M3 in this transition competes 
Mi and E 2 . 

The other phenomenological study [28| applying the 
nonrelativistic QCD (NRQCD) gives the transition width 
as 



8 "l fc | 3 / 2 , 2 , 2n /ooA 

675^ (ai+a2 + fl3) (33) 
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where eti, 02, and 03 are equivalent to the standard form- 
factors Mi, E 2 , and M 3 up to a constant factor and are 
calculated explicitly in NRQCD. By comparing this equa- 
tion with Eq. ([14"|). the factor is approximately 0.33 GcV, 
say, Fi ~ (0.33 GeV)^. Thus, their work gives the pre- 
dictions 

Mi - 0.026 - 0.045 GeV, E 2 ~ M 3 ~ -0.13 GeV, (34) 
which are also in reasonable agreement with our results. 

IV. CONCLUSION 

We calculate in the quenched approximation the mass 
of J PC = 2 ^ charmonium rj C 2, as well as its radia- 
tive transition width to J/ip. The computations are car- 
ried out on two anisotropic lattices with different lat- 
tice spacings, such that the lattice artifacts can be con- 
trolled to some extent. As a calibration, we calculate 
first the spectrum of the lowest-lying charmonia, such as 
IS and IP states, and reproduce the physical pattern of 
the spectrum. In addition, we calculate the transition 
width of Xc2 — > "fJ/4> and get the result 361 ± 9keV, 
which is in good agreement with the experimental value 
of 380 ± 30 keV. Both of these facts manifest the small 
systematic uncertainties due to the quenched approxima- 
tion and the finite lattice spacings. 

There are two states observed in the 2 h channel, 
with masses 3.80(3) GeV and 4.43(8) GeV. The lower 
state has a mass similar to that of the well-established 
■0(3770), which is always assigned to be mainly the l 3 Di 
charmonium, and therefore can be naturally identified 
as the conventional 1 1 Z?2 charmonium 77^. Obviously, it 
is about 70 MeV lower in mass than X(3872), and this 
difference cannot be attributed to the systematic uncer- 
tainties of our work. 

As for the transition rate of rj C 2 — > J/ip'J, we get a small 
partial width of roughly 3.8(9) keV, which is in agree- 
ment with the phenomenological studies. Taking the 
branch ratio Br (X (3872) -> J/ipj) > 0.9% (BABAR) 
or 0.6% (Belle), the full width of V(3872) is estimated 
to be < 420 — 630 keV, which is smaller than, but 
not in contradiction with the experimental upper limit 
Tx < 1.2 MeV. Obviously, a reliable calculation of the 
partial width r\ C 2 — > ip'^f is also crucial for the 1 7?2 char- 
monium assignment of X(3872), but unfortunately there 
are difficulties in the unambiguous extraction of excited 
states on the lattice. However, we can still infer some use- 
ful information from the calculation of r\ C 2 — > J ftp 7- In 
the potential quark model, it is known that 1 Z?2 — > jV is 
a hindered transition with Mi(0) = 0; therefore, the ob- 
served nonzero Mi(0) and the appearance of the higher 



multipoles E2 and M3 can be understood as the rela- 
tivistic correction and the recoil effects of the emission 
of a hard photon (E 1 ~ 0.65 GeV for the final J/tp and 
E 1 ~ 0.11 GeV for ip 1 ). Intuitively, this kind of effect 
for the final tp' can be similar to that for the final J/ip, 
or even milder; thus, the width of r\ C 2 — > 'yip' will be 
suppressed by a kinematic factor of (0.65/0. II) 3 ~ 200 
when compared with the transition rj C 2 — > 'yJ/tp. With 
this fact in mind, the r\ C 2 assignment of X(3872) can 
be ruled out if BABAR's observation of Br(X(3872) -> 
7 ^')/Br(A(3872) -> jj/tp) = 3.4 ± 1.4 is confirmed. 
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Appendix A: Multipoles decomposition of Xc2 <-> J/ip 

For the convenience of readers, the details of the multi- 
pole decomposition of matrix elements of the electromag- 
netic current j At (0) between a vector V state and a 
2 ++ tensor state are described here following Ref. [35| . 
The most general Lorentz-covariant decomposition with 
P and C parity invariance is 

(V(p v ,Xy)\r(0)\T(p T ,X T )) = a(Q 2 )A^ + b(Q 2 )B^ 
+ c(Q 2 )C" i + d T {Q 2 )D» + d v (Q 2 )D» 
+ MQ 2 )F$ + fv(Q 2 )F$ (Al) 

with the definitions 

A? = e»»(p T ,\TK(p v ,\ v ); 

B» = e^(p T ,\ T )pV( e *(p v ,\ v )-p T ); 

C» ee e^{p v ,X v )(e^(pr,X T )plpJ); 

Dt = P^(PT,XT)e* a (p v ,Xv)p^); 

Dy = P^(e^(pT,X T )e* a (pv,Xv)pJ); 

F-t = P^{pt,Xt)pIp V W{pv,Xv)-Pt); 

Fy ee p^(p T7 X T )p^)(e*(p v ,X v ) ■ p T ). 

(A2) 

On the other hand, the matrix elements (V|j^|T) can be 
also expressed in terms of the helicity amplitudes: 
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(V\r\T)e;(X^±) = 



'2fcH 


- 1 


2J- 


- 1 


'2fcH 


- 1 


2J- 


. 

- 1 



(V|j"|T) e *(A 7 = 0) = E \/ |jTT Cfe 5(l + (-l) k SP){kO; J'X\JX), (A4) 



where Ek, M^, and Ck are multipole amplitudes, and SP 
is the product of the P parity of the initial (T) and finial 
(V) states, taking the value SP = — 1 for the 2 + — » 1~ 
transition. These are actually the helicity selection rules. 
An additional constraint comes from the conservation of 
the vector current, 



(v\r\T) q , = d,(v\r\T) = o. 



(A5) 



With the constraints of Eqs. (|A31 IA4| IA5|) . we can solve 
a(Q 2 ), b(Q 2 ), c(Q 2 ),... in the rest frame of the initial 
state with the spatial momentum of the photon parallel 
to the z axis. [The polarization vector of the photon takes 
(1, 0, 0, 1).] Thus, we can get the expressions in terms of 
Ek(q 2 ), Ck{Q 2 )- After that, the general expression of the 
form factor a, b, c, . . . can be obtained by carrying out a 
general Lorentz transformation. For the case of the 2 + — > 
1" transition here, Eq. (|A3|) provides three independent 
equations with respect to the three different helicities of 
the vector state, and therefore gives the relations between 
[a(Q 2 ),b(Q 2 ),c(Q 2 )} and £i(Q 2 ), M 2 (Q 2 ), E 3 (Q 2 )} as 



with Q = (jot • Pv) 2 ~ Tn^my. The constraints from 
Eqs. (|A4[) and (|A5[) with plus/minus helicity of the vector 
meson can fix the parameters dy and g?t- Furthermore, 
fv and /t can be derived from Eqs. (|A4|I and (|A5j) with 
zero helicity of the vector meson. As such, the multipole 
decomposition of the matrix elements (VIj^IT) can be 
expressed finally as 



(V(p v , \ v )\f(Q 2 )\T(p T , A T )) = a^E^Q 2 ) 
+a%M 2 (Q 2 ) + a%E 3 (Q 2 ) + <d(Q 2 ) + a^C 3 (Q 2 ) 

(A7) 



\/l5 

(zV5Ei(m T m, v -pr-Pv) 

+5M 2 (m I m v + Pt-Pv) 
- VZE 3 {Am T m v +p T .p v ) 

2 V5E 3 + 4M 2 
m T = , 

2-y/3ft 



(A6) with the functions af defined by 
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a, 



.4" 



2 



1 

15 



'If 



— 2rny)Fj! + mj(tJ — mvjfyj 
o 

(|cD + 4m y )F^) 



■6^ + -— - m T D v 



■ ( - (Co 2 + Auimy + \rriy)F% + m,^ 



V? = 4 - 



3 TOj 1 



2 TOy 



(my- — LurriT)Dj, + {m T — Cbmr)D v pFC^ — (i m v — wror)^ + (ntr — £jmT)Fy 

(m v — wmr)Dy + (to^ — urnir)Dy ^""(^ I mi/ ) (( TO v — u>rnx)F^ + (m^ — ojm^Fy^J . 



(A8) 



where w = and w = pt-Pv i myiriT- 

Appendix B: Multipoles of ?7 C 2 -o J/^ 

The case of r; C 2 — > J/V' is sightly different from that 
above. The most general Lorentz-covariant decomposi- 
tion with P and C parity invariance is 

(V(pv,XvW(0)\ V 2(pt,Xt)} = a{Q 2 W + HQ 2 )B» 
+ c{Q 2 )C» + d{Q 2 )D» + e(Q 2 )E» (Bl) 

with 

A" = e^et(fv,X v )pJ4(PrAT)pJ, 

B» = e^pJet(p v ,X v )p^e^p T ,X T ), 

c , = e ^pT p v e P^ T] X T )e* (py,X v ), 

D" = e^plp T p e*(p v ,X v )e^(p r ,X T )p^ > 

E , = e ^ p y p T e ^ Tj x T )p^* a (p v ,Xv)p T a . (B2) 



In fact, there exist another three Lorentz-covariant struc- 
tures At, Ey, and Et- 

E v = e a ^p T a e}{p v ,X v )p v p e^{p T ,X T )pJpy, 

A T = e^et(py,X v )p T p e^PT,X T )p^, 

E T = e a ^ple}(p v ,X v )p^e^(p T ,X T )pJp T , (B3) 
which, however, arc not independent and can be ex- 
pressed in terms of the functions in Eq. (|L32|) : 

A T = -B^-C 

E T = m 2 T A t ' l +p v .p T B t " +p v .p T C t " +E" 

E v = py PT A^ + m 2 v B^ + m 2 yC^ + D^. (B4) 

So they do not appear in the decomposition. Based on 
this, one can follow the similar procedure of the case of 
Xc2 J/ip to derive the related multipole decomposi- 
tion. The constraints of decomposition are similar to 
Eqs. (|A3IA5I) while SP = 1. Finally, one can get the 
result 
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(V(p v> X v )\r(0)\ m (pT,X T )) = 
' 1 iQ ~ 1 [ - Vl5C» + ^VT5E»(-m v m T + Pt-Pv) 



* E 2 (Q 2 ) 



3nVa 
i M 3 (Q 2 ) 



(mym T +pt-Pv) 



* c 2 (0 2 ) 



- 2v / 15C A1 + - [5VT5D tl m^ + 2VT5E t *(4m v m T + Pt-Pv) 
A^rriT + B^rriT + C^mx + tz ^mr^ - Pt-Pv) + E^m T {m\ - pt-Pv) 



(B5) 



Appendix C: Form factor as function of Q 2 or Q 

Transition form factors are always expressed as Lorentz 
scalar functions of the squared momentum transfer, Q 2 = 
— g 2 = — {q% — g/) 2 , where pi and pf refer to the four- 
momenta of the initial and final particles, respectively. 
However, if one looks into the Lorentz decomposition 
[Eq. (|Aip ] and the multipole decomposition [Eq. (|A7I) ] 



for the Xc2 — > J/fp transition matrix elements, one 
can find that the quantity fl = (pi.pf) 2 — m 2 m 2 , with 
Pi-Pf — (mf + m 2 + Q 2 )/2, is also an interesting Lorentz- 
invariant kinematic variable. According to Eq. (|A6|) , the 
multipole amplitudes E\, M2, and £3 can be reversely 
expressed in terms of the form factors a, 6, c as, 



M 2 (0) 



6c my — 3a tot 



-a 

3 4Vl5mT 
2c my — a tot + 46 m v mT 

—2c my + (a + 2& to 2 /) Tot 



15tot 




(Cl) 



Obviously, these are polynomials of the variable 
f2/ (to 2 / to 2 ") with the coefficients the combinations of the 
form factors a, 6, and c. The physical meaning of the 
above expression can be understood in the rest frame of 
the decaying particle (T here), where Q/(m v m T ) = v 2 
with v = [pV|/TOy ~ 0.16 for the \c2 —tJ/ip transi- 
tion. It is clearly seen from these expressions that the 
Ei transition is dominant, while Mi and E3 transitions 
are suppressed by a factor of v 2 ~ 0.026. On the other 



hand, M 2 (Q 2 ) and E 3 (Q 2 ) should be zero at v = 0, or 
equivalently, Q 2 = — (my — my) 2 , as is confirmed by our 
simulation results (seen in Fig. @|. It is found in the cal- 
culation that M2 and E3 are consistent with zero when 
both the initial and final state are at rest. 

As for the rj C 2 — > J/ip transition, the multipole ampli- 
tudes Mi , E2 and M3 can similarly be expressed as poly- 
nomials in Q/(TOym r ), where the coefficients are also the 
combinations of the form factors in Eq. (|B1|) : 
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rriT 



M 3 (n) 



— (a my + a tot — 2c tot) 



-2a mv + 3c tot 



16V15 



tijn / / b~y 



2a TOy — 3c tot — 4c? to^tot + 6e mym T 



SI 



4V15 



n 



fciyi / / b~y 



£ 2 (Q) = 

rriT 



3 , s 2a TOy — c tot — 4d to^tot + 2e myTO T 



— (a Tot — a wiy) + 
4 



-2a my + c tot 



16%/3 



SI 



brj-\ lib y 



o 



brj-\ lib -^y 



4^3 

3" 



rriT 
+ 01 



-a my — c Tot + 2d TOy tot + 2e mym 



CI 



by ' ' v'J^ 



n 



brj-i I I by 



n 



by ' i 



a my + c tot 



4V15 



n 



(C2) 



In the rest frame of the decaying particle r] C 2 (denoted physical implication of these expressions has been dis- 
by T here), the prefactor \/U/m T = \p v \ is exactly the cussed in context and is omitted here, 
requirement of the P-wave decay of i] C 2 — > jJ/ip- The 
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